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Abstract 

We consider Sakai-Sugimoto like models consisting of Dq-Dp-Dp-branes where iVf flavor Dp and 
Dp-branes transversely intersect N c color Dg-branes along two (r + l)-dimensional subspaces. For 
some values of p and q, the theory of intersections dynamically breaks non-Abelian chiral symmetry 
which is holographically realized as a smooth connection of the flavor branes at some point in the 
bulk of the geometry created by N c Dg-branes. We analyze the system at finite temperature and 
map out different phases of the theory representing chiral symmetry breaking and restoration. 
For q < 4 we find that, unlike the zero-temperature case, there exist two branches of smoothly- 
connected solutions for the flavor branes, one getting very close to the horizon of the background 
and the other staying farther away from it. At low temperatures, the solution which stays farther 
away from the horizon determines the vacuum. For background D5 and D6-branes we find that 
the flavor branes, like the zero temperature case, show subtle behavior whose dual gauge theory 
interpretation is not clear. We conclude with some comments on how chiral phase transition in 
these models can be seen from their open string tachyon dynamics. 



1 Introduction, summary and conclusions 



A very interesting holographic model of QCD which realizes dynamical breaking of non-Abelian 
chiral symmetry in a nice geometrical way is the Sakai-Sugimoto model [1]. In this model, one starts 
with N c D4-branes extended in (x°a; 1 x 2 a; 3 a; 4 )-directions with x A being a circle of radius R. The low 
energy theory on the branes is a (4 + l)-dimensional SU (N c ) SYM with sixteen supercharges. To 
break supersymmetry, anti-periodic boundary condition for fermions around the x 4 -circle must be 
chosen. To this system, Nf D8 and iVf D8-branes are added such that they intersect the D4-branes 
at two (3 + l)-dimensional subspaces R 3, , and are separated in the compact :r 4 -direction by a 
coordinate distance of £q = irR. In other words, the D8 and D8-branes are located asymptotically 
at the antipodal points on the circle. The massless degrees of freedom of the system are the 
gauge bosons coming from 4 — 4 strings and chiral fermions coming from 4-8 strings. (Before 
compactifying x 4 , there are also massless adjoint scalars and fermions coming from the 4 — 4 
strings. These modes become massive upon compactifying x 4 and choosing anti-periodic boundary 
condition for the adjoint fermions; fermions get masses at tree level whereas scalars get masses due 
to loop effects.) The fermions localized at the intersection of D4 and D8-branes have (by definition) 
left-handed chirality and the fermions at the intersection of D4 and D8-branes are right-handed. 
The U(N{) x U(N{) gauge symmetry of the D8 and D8-branes in this model is interpreted as the 
chiral symmetry of the fermions living at the intersections. At weak effective four-dimensional 't 
Hooft coupling A4 the low energy theory contains QCD but this is not the limit amenable to analysis 
by the gauge-gravity duality. At strong-coupling (A4 3> 1), however, the theory is not QCD but 
can be analyzed using the gauge-gravity duality [21 El HI [5] , and it has been suggested that it is 
in the same universality class as QCD. In lack of any rigorous proof for this universality, the best 
one can do is to check whether this model at large N c and large four-dimensional 't Hooft coupling 
exhibits the key features of QCD; namely, confinement and spontaneous chiral symmetry breaking. 
In fact, it apparently does pQ. By considering Nf D8 and iVf D8-branes as probe "flavor" branes 
in the near-horizon geometry of N c "color" D4-branes and analyzing the Dirac-Born-Infeld (DBI) 
action of the flavor branes in the background, one observes that at some radial point in the bulk the 
preferred configuration of the flavor branes is that of smoothly-connected D8 and D8-branes. This 
geometrical picture is interpreted as dynamical breaking of U(N{) x U(N{) chiral symmetry (where 
the branes are asymptotically separated) down to a single U(Nf) (where the branes connect). The 
model also shows confinement [H [7] . 

Although choosing £q = ttR makes calculations a bit simpler, there is no particular reason to 
consider the flavor branes to be asymptotically located at the antipodal points of x 4 . In fact, the 
£q <C ttR limit of the Sakai-Sugimoto model is very interesting in its own right. By analyzing the 
Iq <C ttR limit, or equivalently the R —* 00 limit, it was realized in [8] that the theory at the 
intersections can be analyzed both at weak and strong effective four-dimensional 't Hooft coupling 
A4. At weak-coupling the model can be analyzed using field theoretic methods and in fact is a non- 
local version of Nambu-Jona-Lasinio (NJL) model [9]. At strong-coupling it can be analyzed by 
studying the DBI action of the flavor branes in the near-horizon geometry of the color branes and 
exhibits chiral symmetry breaking via a smooth fusion of the flavor branes at some (radial) point in 
the bulk. A nice feature of this model is that the scale of chiral symmetry breaking is different from 
that of confinement [8], and one can completely turn off confinement by taking the R — > 00 limit. 
Therefore the £q <C ttR limit of the Sakai-Sugimoto model provides a clean holographic model of 
just chiral symmetry breaking without complications due to confinement. 
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The finite-temperature analysis (at large N c and large 't Hooft coupling A4) of the Sakai- 
Sugimoto model as well as the holographic N JL model [8] , was carried out in |10l [TT] . Putting the 
flavor branes as probes (iVf <C N c ) in the near-horizon geometry of N c non-extremal D4-branes and 
analyzing the DBI action of the flavor branes, one obtains [101 [TT] that at low temperatures (com- 
pared to Iq 1 ) the energetically-favorable solution is that of smoothly-connected D8 and D8-branes 
which, like its zero-temperature counterpart, is a realization of chiral symmetry breaking. At high 
enough temperatures, on the other hand, the preferred (in the path integral sense) configuration is 
that of disjoint D8 and D8-branes, hence chiral symmetry is restored. Also, when is compact 
and £0 < ttR, there exists an intermediate phase where the dual gauge theory is deconfined while 
chiral symmetry is broken |10j . 

It is certainly interesting to explore whether the holographic realization of chiral symmetry 
breaking and restoration is specific to a particular model such as the aforementioned ones, or generic 
in the sense that other intersecting brane models will realize it, too. To search for genericness (or 
non-genericness) of chiral symmetry breaking in intersecting brane models, a system of Dq-Dp- 
Dp-branes was considered at zero temperature in |12] where the color Dg-branes are stretched 
in non-compact (x°x 1 . . . x 9 )-directions. The flavor Dp and Dp-branes intersect the color branes 
at two (r + l)-dimensional subspaces K r ' . Without flavor branes, the low energy theory on the 
color branes and whether it can be decoupled from gravity (or other non-field theoretic degrees 
of freedom) with an appropriate scaling limit was analyzed in [13]. The low energy theory on the 
Dg-branes is asymptotically free for q < 3, conformal for q = 3, and infrared free for q > 3. While 
for q < 5, there always exists a scaling limit for which the open string modes can be decoupled from 
the closed string modes, there exists no such limit for q = 6. With the flavor branes present, an 
analysis was carried out in [12] where the behavior of the model for both small and large effective 
't Hooft coupling A g- ~ \ q +i(.Q~ 9 (£q is the asymptotic coordinate distance between Dp and Dp- 
branes, and A q +i is the 't Hooft coupling of the (q + l)-dimensional theory on the Dg-branes.) was 
considered. Taking the N c — > 00 limit while keeping A e fr fixed and large, amounts, in the probe 
approximation, to putting the flavor branes in the near-horizon geometry of N c color branes. (See 
[121 [T5] for the validity of the supergravity analysis in these models.) Determining the shape of 
the flavor branes (relevant for chiral symmetry breaking) by analyzing their DBI action in the 
background geometry, it was observed that [12] for q < 4 there always exists a smoothly-connected 
brane solution which is energetically favorable. For a subclass of these general intersecting brane 
models, namely those for which q + p — r = 9, the above-mentioned connected solutions can be 
identified with the U(N{) x U(N{) chiral symmetry being spontaneously broken. Following [12] 
we will call the brane models for which q + p — r = 9 as transversely-intersecting brane models 
and their intersections as transverse intersections. For q = 5, there exists no connected solution 
except when £q takes a particular value. For this particular value of £q which is around the scale 
of non-locality of the low energy theory on the D5-branes (which is a little string theory; see [14] 
for a review of little string theories), there is a continuum of connected solutions whose turning 
points can be anywhere in the radial coordinate of the bulk geometry. All such solutions are equally 
energetically favorable. For q = 6, there is a connected solution but is not the preferred one. The 
more energetically favorable solution, in this case, is that of disjoint branes. Due to the lacking of 
an appropriate decoupling limit for q = 6, it is not clear whether one can realize such a solution as 
a phase for which chiral symmetry is unbroken. 



3 



Summary and conclusions The purpose of this paper is to investigate some aspects of non- 
compact transversely-intersecting D-brane models at finite temperature, in particular the number 
of solutions, their behavior, and whether or not such solutions can be identified with a chirally 
broken (or restored) phase of the dual gauge theory living at the intersections. The main reason 
for us to consider transverse intersections is that in the probe approximation the generalization 
of the Abelian U(l) x U (1) chiral symmetry to the non-Abelian case is straightforward: one just 
replaces Nf = 1 with general iVf, and multiplies the flavor DBI action by Nf. This is not the 
case in other holographic models of chiral symmetry breaking and restoration^. Note that the 
Sakai-Sugimoto model p], its non-compact version [8], and the D2-D8-D8 model analyzed in [T7] 
which holographically realizes a non-local version of the Gross-Neveu model |18| are examples of 
transverse intersections. 

This paper is organized as follows. In section two we first review the general set up of the 
transversely-intersecting Dg-Dp-Dp-branes, then consider the system at finite temperature at large 
N c and large effective 't Hooft coupling A g. There are two saddle point contributions (thermal and 
black brane) to the bulk Euclidean path integral, and each can potentially be used as a background 
geometry dual to the color theory. Comparing the free energies of the two saddle points, we 
determine which one is dominant (has lower free energy) for various g's. For q ^ 5, the dominant 
saddle point is the black brane geometry, whereas for q = 5, the thermal geometry is typically 
dominant. In section three we present the solutions to the equation of motion for the DBI action 
of the flavor branes placed in the near horizon geometry of black Dg-branes. By a combination 
of analytical and numerical techniques, we show that, unlike the zero-temperature case, for £q/(5 
less than a critical value, there exist generically two branches of smoothly-connected solutions (and 
of course, a solution with disjoint branes) for q < 4. (/? is the circumference of the asymptotic 
Euclidean time circle, and is equal to the inverse of the dual gauge theory temperature T.) Note 
that some of these branches were previously missed in the literature. One branch which we will call 
the "long" connected solution gets very close to the horizon of the black Dq-branes whereas the 
other one named "short" connected solution stays farther away from it. Beyond the critical value, 
the flavor branes are "screened" and cannot exist as a connected solution. The situation is, however, 
different for q = 5,6. For q = 5 and T < (27ri?5 + i)~ 1 , where -R5+1 denotes the characteristic radius 
of the D5-brane geometry, the flavor branes must be placed in the near horizon geometry of thermal 
D5-branes. Like the zero temperature case, we find that there exists an infinite number of connected 
solutions when £q is around the non-locality scale of the low energy theory on the D5-branes. This 
scale is set by the (inverse) hagedorn temperature of D5-brane little string theory. There are no 
connected solution for other values of £q, though. For T = (27ri?5 + i) _1 the flavor branes should be 
considered in the near horizon geometry of black D5-branes. In this case there always exists one 
connected solution, as well as a disjoint solution, for small £0 (compared to (27ri?5 + i) _1 ). In fact, 
at this particular temperature there is one connected solution for ^o' s much less than -R5+1. The 
number of solutions for £q beyond the non- locality scale (~ R5+1) depends on the dimension of the 
intersections. For four-dimensional intersections, there are two connected solutions up to a critical 
value whereas for two-dimensional intersections there is no connected solution. Of course, for 
both two- and four-dimensional intersection there is always a solution representing disjoint branes. 
Having determined the flavor brane solutions in the background of thermal and black D5-branes, 
it is not clear whether or how these solutions represent a chirally-broken or restored phase in the 

For non-transverse intersections, namely q + p — r 7^ 9, one can identify a symmetry in the common transverse 
direction as a chiral symmetry for the fermions of the intersections |15] . In some cases the generalization to low-rank 
non-Abelian chiral symmetry is possible [16], but such generalizations are not generic. 
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dual theory, mainly because in the geometry of D5 — branes, there are modes (non-field theoretic) 
which cannot totally be decoupled from the dual field theory degrees of freedom. Lastly, for q = 6, 
independent of what value £o/(3 takes, there always exist one connected solution (and a disjoint 
solution). 

In section four, we map out different phases of the dual gauge theories and determine whether 
or not there is a chiral symmetry breaking-restoration phase transition. We do this by comparing 
the regularized free energies of the various branches of the solutions found in section three. For 
q < 4 the short solution is preferred to both long and disjoint solutions at small enough tempera- 
tures (compared to £q 1 ), hence chiral symmetry is broken. At high temperatures, however, chiral 
symmetry gets restored and this phase transition is first order. For q = 5 and T < (27ri?5 + i) _1 , 
the infinite number of connected solutions are all equally energetically favorable, and each one of 
them is preferred over the disjoint solution. For T = (27ri?5+i) -1 , we find that for small enough 
£o/(2vri?5 + i) the disjoint solution is preferred. For larger values of £q, there is no connected so- 
lution so the disjoint solution is the vacuum. As we alluded to earlier, it is not clear to us that 
the preferred solutions of the flavor branes in the background of color D5-branes can be associated 
with different phases of the dual theory. For q = 6, although we find that the disjoint solution is 
always preferred and there is no phase transition, there is no clear way to associate this solution 
with unbroken chiral symmetry in the dual field theory. This is because the is no decoupling limit 
that one can take to separate the gravitational degrees of freedom of those of the dual theory. 

Section five is devoted to a brief analysis of the number of solutions and their energies for 
transversely- intersecting Dq-Dp-Dp-branes with compact x q . In section six we speculate how the 
order parameter for chiral symmetry breaking can be realized in finite-temperature transversely- 
intersecting D-branes by including the thermal dynamics of an open string tachyon stretched be- 
tween the flavor branes, and how it may depend on temperature. Finally, in the appendix we 
present detailed calculations for the free energies of the near horizon geometries of color Dg-branes 
(with the topology of either S 1 x S 1 or S 1 x R in the t — x q submanifold) to determine the dominant 
background geometry (either thermal or black brane) at low and high temperatures. 



2 Transverse intersections at finite temperature 



We start this section by reviewing first the general setup for transverse intersections of Dq-Dp-Dp- 
branes and identifying the massless degrees of freedom at intersections. We consider the system at 
finite temperature in the large N c and large 't Hooft coupling limits. Since there is more than one 
background, we determine the one with the lowest free energy and consider that as the background 
dual to the color sector of the dual theory at finite temperature. We then write the equation of 
motion for the flavor branes. This section is followed in the next section by an analysis of the 
solutions to the equation of motion as a function of the dimensions of the intersections r, as well 
as q. (For transverse intersections p, the spacial dimension of the flavor branes, is determined once 
r and q are given.) 
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2.1 General setup 



Consider a system of intersecting Dg-Dp-Dp-branes in flat non-compact ten-dimensional Minkowski 
space where N c Dg-branes are stretched in . . . x (? )-directions, and each stack of Nf Dp and 

Nf Dp-branes are extended in (x°x 1 . . . x T )- and (x q+l . . . x 9 )-directions. The Dp- and Dp-branes 
are separated in the a^-direction by a coordinate distance £q and intersect the Dq-branes at two 
(r + l)-dimensional intersections 



Dq x 
Dp x 
Dp x 

Using T-duality one can determine the massless degrees of freedom localized at the intersection. 
It turns out that for transverse intersections, q + p — r = 9, the massless modes which come from 
the Ramond sector in the p — q strings are Weyl fermions. These fermions are in the fundamentals 
of U(N C ) and U(Nf). The massless modes at the other intersection are also Weyl fermions which 
transform in the fundamentals of U(N C ) and U(Nf) of the Dp-branes. 



x 
x 
x 

X 



X' 

X 
X 
X 



X 
X. 



(1) 



One way to put the above system at finite temperature (at large iV c , large effective 't Hooft 
coupling A e ff, and in the probe approximation Nf <C N c ) is to start with the geometry of black 
Dg-branes as background. The Euclidean metric for this geometry is 

ds 2 = ( -^-) ~ ( f(u)dt 2 + dx 2 ) + ( ~~ (£L + u 2 dnl_ q ) , (2) 

with 



/w=i- (^r, (3) 

where in the metric dttg_ q is the line element of a (8 — g)-sphere with a radius equal to unity, and 
Rq+i, which denotes the characteristic radius of the geometry, is given by 

K+l = (2^f- q T {^) 9s Nfr = ^(V^) 9 -^r^y q+1 NJl^, (4) 

where g s is the string coupling. The Euclidean time is periodically identified: t ~ t + /3, where j3 
is equal to the inverse of the temperature T of the black branes. In ([3]) the horizon radius ut is 
related to (3 as 



47r \R q+ iJ ut 

The relationship between ut and (3 comes about in order to avoid a conical singularity in the metric 
at u = ut- Note that for black D5-branes (3 is independent of ut, and equals 27ri?5+i. 

Also, the dilaton (j> and the g-form RR-flux F q are given by 

/ u \ L ^- 3 ^ 2vriV c 

e=9s {R^) ' F «=V^ e8 -- (6) 
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where Vs- q and es- q are the volume and the volume form of the unit (8 — g)-sphere, respectively. 



There is, however, another background with the same asymptotics as ([2]) which may potentially 
compete with the aforementioned background. The metric for this (thermal) geometry is 

7-q T-q 

d8*=(j^£) 2 (d^ + d ^ + (^_^ 2 ^ +u 2 dSl 2_^ f (7) 

with the Euclidean time t being periodically identified with a period j3 = T _1 . Unlike the black 
brane geometries ([2]) , (3 could take arbitrary values in the thermal geometries ([7]) . The dilaton and 
the g-form RR-flux are the same as ([6]). 

In the appendix we have calculated the free energies of both thermal and black brane geometries. 
Except for q = 5, the difference in free energies AS* of the two geometries subject to the same 
asymptotics is given by 

A.S = Sthcrmal — Sblack brane = 2 ^9^T ^ 1 (^) 

9s 

where Vg is the volume of space transverse to the radial coordinate u: Vg = Vol(S 8_9 )Vol(M 9 )Vol(Sjg). 
The volume Vg is measured in string unit l s where for simplicity we set l s = 1. The difference in 
free energies (J8j) shows that the thermal background is less energetically favorable compared to 
the the black brane background (|2]). Thus, there is no Hawking-Page type transition between the 
two geometries which holographically indicates that there is no confinement-deconfinement phase 
transition in the dual theory. The situation is different for q = 5. Semi-classically, once the 
characteristic radius R5+1 is given, the black D5-brane geometry will have a fixed temperature 
T = (2"7ri?5 + i)~ 1 . At this specific temperature, there are two saddle points contributing to the 
(type IIB) supergravity path integral: thermal and black D5-brane geometries. The difference in 
free energies of the two saddle points is (see the appendix for more details) 

5'thermal D5 — 'S'black D5-brane = — T-Vol(S 3 )Vol(M 5 )i25 + itty, at (3 = 2irR 5+ i, (9) 

9s 

showing that the black brane geometry is the saddle point with lower free energy. However, at 
temperatures other than (27ri?5_|_i ) — 1 , the thermal geometry is the only saddle point, although due 
to the hagedorn temperature of the D5-brane theory, one should only consider temperatures less 
than (2vri? 5+ i)- 1 . Thus, for T < (2ttR 5+ i) 1 we use the thermal D5-brane geometry as background. 



2.2 Flavor Dp-Dp-branes in black Dg-brane geometries 

We are interested in the dynamics of the flavor Dp and Dp-branes in the background of the black 
Dg-brane geometries. As we alluded to earlier, for q = 5 the thermal geometry of N c D5-branes 
(once the limit of the near horizon geometry is taken) is the background that one should use for 
the dual finite temperature field theory for T < (27ri?5 + i) -1 . We also analyze the dynamics of the 
flavor branes in the background of the black D5-branes in which case it is understood that the dual 
theory is at a fixed temperature T = (27ri?5+i) — 1 . We are interested in the static shape of the 
flavor branes as a function of the radial coordinate u. Therefore we choose the embedding 

t = a°, x 1 = a 1 , ... x r = a r , x q = a^ . . 

u = u{a q ), x q+2 = a q+ \ ... x 8 = a p -\ x 9 = a p , [ } 
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subject to the boundary condition 



u(± — 
v 2 ' 



oo, 



(11) 



where {a ,-- - ,cr p } are the worldvolume coordinates of the flavor branes. This boundary condi- 
tion simply states that the asymptotic coordinate distance between the Dp and Dp-branes is £q. 
Ultimately the stability of such an assumption lies in the large N c limit. 

From now on, we set Nf = 1. As it becomes apparent in what follows, the generalization to 
iVf <C N c is straightforward. The dynamics of a Dp-brane (and a Dp-brane) is determined by its 
DBI plus Chern-Simons action. Solving the equations of motion for the gauge field, one can safely 
set the gauge field equal to zero and just work with the DBI part of the action. After all, it is this 
part of the full action which is relevant for our purpose of determining the shape of the Dp and 
Dp-branes. Therefore, with gauge field(s) set equal to zero, the dynamics is captured by the DBI 
action 



Sdbi = P P J d p+1 a e-^y/det^) 



(12) 



where \i v is a constant, and g a b = Gmn9o,x ' d\,x N is the induced metric on the worldvolume of the 
Dp-brane. For a Dp-brane forming a curve u = u(x q ), the DBI action (|12j) reads 



S BBl = pC(q,r) / d r x dx q f(u) + 



R. 



U \2<5 ,2 
U 



9+1 



(13) 



where 



C(q,r) 



Vol(S 8 -«)i? 3+1 4 



i(,-7)(r-3) 



7 = 2 + -(7-5)(r + l), 
$ = i(ff-7), 



(14) 



and u = du/dx q . For a Dp-brane forming a curve u(x q ) in the thermal D5-brane geometry the DBI 
action is obtained by setting f(u) = 1 in (Tl3l) . The integrand in (Tl3j) does not explicitly depend on 
x q , therefore C — u dC/du must be conserved (with respect to x q ). A first integral of the equation 
of motion is then obtained 



uif(u) f(u) + 



R, 



U \ 2S il 

u 



9+1 



t ' 



(15) 



where uq parametrizes the solutions. We now analyze the solutions of (|15p . 



3 Multiple branches of solutions 



The simplest solution of the equation of motion in (|15p . namely uq = 0, corresponds to x q = 
constant. In order to satisfy the boundary condition (flT|) . one obtains x q = ±Iq/2. So, the uq = 
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solution corresponds to disjoint Dp and Dp-branes descending all the way down to the horizon at 
u = ut- Also, note that the existence of this solution is independent of (3. 

For uo 7^ 0, solving for u yields 

Since the left hand side of (|16p is non-negative, the right hand side of (|16|) must also be non-negative 
resulting in u > maxj-ur,^}, where u* is a possible turning point. Therefore, for allowed solutions 
one must have u > u* > ut- 

The possible turning points are determined by analyzing the zeros of the right hand side of 
([16j). Setting /(«*) = will not result in a valid turning point. So, the other possibilities come 
from solving ii*/(/u*) — Uq = 0, which we will rewrite as follows 

uj - Ku T 25 - ul = 0, (17) 

where 

<r = 7 + !W = 2 + i(7-?)(r-l). (18) 

Note that since (and in fact, for the cases of interest, it is either 1 or 3), a is always a positive 

integer which, combined with the fact that 5 < 0, implies that a < 7. 

We use (fT6|) to relate the integration constant u*, or equivalently uo, to the parameters of 
the theory, namely the (inverse) temperature (3 and the asymptotic distance between the Dp and 
Dp-branes £q. First, rearrange (fTUj) to get 

= R- 5 (u2 - u°u T 25 ) l 2 f (V 25 - u" 25 )^ x (19) 

«/ 1/.* 

(V - u CT n T - 25 - («2 - ulu T 2S )Y^du, 

where in the second line we used (fTT|) to trade uq for u*. Changing to a new (dimensionless) variable 
z = u/ut, CE|) becomes 

x q {z) = - <)3 £ (V 25 - (* 7 - - fe 7 - (20) 

where z* £ (1, 00). Taking the z — > 00 limit, we can relate z* to (3 and ^0 

I = ~ 6 - (zl - £° (z- 25 - - z CT - (z7 - (21) 

As we mentioned earlier, the solutions to the equation of motion are parametrized by possible 
value(s) of the turning point z*. For a fixed £q//3, it is the number of z* which determines the 
number of (connected) solutions. Thus, one has to analyze £q//3 as a function of z* to determine 
the number of solutions for a fixed £0/ (3. 
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3.1 Analytical analysis 



There are regions of z* for which io/(3 as a function of z* can be given analytically. These are the 
z* — ► 1 + and z* ^ 1 regions. For any z*, the integral in (|2ip can be evaluated numerically. The 
numerical results will be presented shortly after the analytical analysis for the two limiting cases 
is given. 



9 = 6 



g = 5 




1 <- z, 



z, » 1 



Figure 1: Behavior of £o//3 versus 2* in two regions of z* — » 1 and z* 3> 1 for various g's. Except for q = 5, these 
plots illustrate the number of connected Dp-branes at low and high temperatures placed in the background of black 
Dg-branes. For q = 5, it is understood that ft is fixed; j3 = 2nR$ + i, and different solutions is obtained by varying £q. 



First consider the z* — ► 1 + limit. Taking z* = 1 + e where 0<f<l, (|2ip is approximated by 
4 5 



— ~ \/-25e I ( 

P 7T J 1+e 



-25 



l)"2( z T- z CT + 2(5e)-2dz. 



(22) 



Ignoring some numerical prefactors, the behavior of (|22|). is approximated by y/e f°° z s zdz for 
large values of z. Since 5 — ^ < —2 for all g and r of interest, one has \/e f°° z 5_ i c?z ~ y^e. On the 
other hand, when z approaches z* such that z — z* > 0, we define z = z* + x with < x <C 1, and 
expand out the integrand of (|22|) around x. We get 



£0 




==Ve / + e)) 2dx 

TT\J — Id Jo 
1 



~ V-2(5 v / elog Ve, 

7T 

indicating that the leading behavior of £o//3 in the z* — > 1 limit is — y^log y^e. 
For the z* 3> 1 region, we can approximate (f2Tj) by (recall 5 < and <r < 7) 

! dz + • • • 
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where • • • represents terms subleading in z* , and in the second line in (|24p we have changed the 
variable from z to y = z/z*. Thus, aside from a numerical factor, the z* — > oo limit of (|2ip reads 

| - (25) 



This expression is identical to the one derived for the zero temperature case in |12| . This resemblance 
is not accidental because the large z* limit corresponds to having a turning point very far away 
from the horizon of the background geometry. The results obtained for this limit should then match 
those derived for the zero temperature case. An interesting feature of (]25h is that for the black 
D5-branes £o/(2ttR^ + i) is independent of z* and approaches a constant value of l/(r + 3). This 
value has been argued in |12j to be around the scale of non-locality of the low energy effective 
theory on D5-branes. The analysis for the dynamics of the flavor branes placed in the thermal D5- 
brane geometry is the same as the analysis when they are placed in the zero temperature D5-brane 
geometry. The zero temperature analysis has already been done in [12j where it was found that 
there exist an infinite number of connected solutions for one specific value of £q = 2-7ri?5+i/(r + 1), 
and none for other £q's. 



Analyzing (|2ip for the two regions of z* — > 1 + and z* 3> 1, the minimum crude conclusion that 
one can draw is that for small enough £o/P there exist two connected solutions (one closer to the 
horizon which we will name "long" connected solution, and the other farther away from it named 
"short" connected solution) for q < 4 and only one curved solution for q = 6. For q = 5 and for 
T < (27ri?5 + i)~ 1 , there exists an infinite number of connected solutions for just £q = 27ri?5 + i/(r + l) 
and none for others. On the other hand, for T = (2-7ri?5 + i) _1 , there is only one connected solution 
given that £q <C 2-7ri?5 + i. The analysis for the two z* regions has been summarized in Figure [TJ 



3.2 Numerical analysis for the number of solutions 



For generic values of £o/P, the allowed number of solutions can be determined by numerically 
integrating OTj) and plotting £o/j3 versus z*. The results for various intersections and Dg-branes 
are shown in Figure 2. Although we have plotted £o/f3 versus z* for z* £ (1,50), the qualitative 
behavior of the plots stays the same if one considered larger values of z*. As we will describe 
below, the number of curved solutions depends on what configuration is being considered. Note 
that regardless of the configuration, there always exists a disjoint solution. In what follows in the 
rest of this subsection, when we say there exist one or two solutions for a particular system we have 
connected solutions in mind. 



(3 + l)-dimensional intersections In this case there are three allowed D-brane configurations 
which are transversely intersecting. These are the D4-D8-D8, D5-D7-D7 and D6-D6-D6 configura- 
tions shown on the top row in Figure [2J 

For the D4-D8-D8 configuration (shown on the upper left corner in Figure 2), there is a critical 
value of (£o/P)cr ~ 0.17 beyond which there exists no connected solution for the flavor branes. 
Below this critical value there are two connected solutions which we earlier called the "short" and 
the "long" connected solutions. The existence of these two types of solutions and the critical value 
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of 0.17 were already noted by the authors of [TT] in the their analysis of the holographic NJL 
model at finite temperature. We will see in the next section that the short solution is always more 
energetically favorable to the long one. Since there also exists a disjoint solution, determining the 
chirally-broken or chirally-symmetric phase of the dual field theory (which is a non-local version 
of the NJL model [8]) is just a matter of comparing the free energies of the disjoint and short 
connected solution. This will be done in the next section. 

The analysis for the D5-D7-D7 and D6-D6-D6 cases are, however, more subtle, and the holo- 
graphic interpretation of the solutions is less transparent for reasons to be mentioned below. For 
the D5-D7-D7 configuration at temperature T = (27ri?5 + i) _1 , there are two critical values of 
(£o/(2vri?5 + i)) cr 0.168 and 0.175. For ^o/(2vri?5+i) < 0.168 there is always one connected solu- 
tion, for 0.168 < £q/(2ttR 5+ i) < 0.175 there are two, and for £ /(2ttR 5+ i) > 0.175 there exists none. 
Comparing these results to the ones obtained for the zero-temperature D5-D7-D7 configuration, 
one observes that while at zero temperature |12| . there are either an infinite number of solutions 
or none, at T = (27ri?5 + i) -1 there are different numbers of solutions (one, two, or none) depend- 
ing on what values £q take; see Figures 2 and 3. For D5-D7-D7 configuration at temperatures 
T < (2-7ri?5_|_i) _1 , the situation is the same as the zero temperature case. There exist connected 
solutions only for a specific value of £o = ttR^ + i/3. Indeed, there are an infinite number of such 
solutions; see Figure 3. 

For D6-D6-D6 configuration, the situation is simpler. For any £q/ (3 there always exists one and 
only one connected solution. As we will see in the next section all these connected solutions are less 
energetically favorable compared to disjoint solutions. The observation that for color D6-branes 
at finite temperature there is always one connected solution, hence no "screening length", is in 
accord with the fact that the near horizon geometry of D6-branes cannot be decoupled from the 
gravitational modes of the bulk geometry [13j . 



(1 + l)-dimensional intersections In this case, there are five allowed configurations (shown on 
the second and third rows in Figure 2), namely D2-D8-D8, D3-D7-D7, D4-D6-D6, D5-D5-D5 and 
D6-D4-D4. Configurations with q < 4 show similar behavior. For q < 4 there always exists a 
critical (£o/(3) C r beyond which connected solutions cease to exist. The critical value obtained from 
Figure 2 is (£ /(3) CI « 0.225,0.223,0.227 for color D2, D3 and D4-branes, respectively. Below these 
critical values there are two solutions, the short and long solutions^. 

For the D5-D5-D5 system at T = (27ri?5 + i) _1 , below a critical value of (4)/(27ri?5 + i)) cr 0.251 
there is always one solution whereas above this value connected solutions do not exist. Note 
the difference (depicted in Figure 3) of this case with the D5-D7-D7 system: for the D5-D5-D5 
configuration, there is no range of £q/ {2-kR^i) for which there exist the short and long solutions. 
Like the D5-D7-D7 configuration, for the D5-D5-D5 system at T < (27ri?5 + i)~ 1 , there is an infinite 
number of connected solutions when 2£q = ttR^ + i and none for other £o's. 

For the D6-D4-D4 configuration, there is always one connected solution for arbitrary values of 
£o//3- Notice that the behavior of the flavor branes in the geometry of color D5 and D6-branes with 
(l + l)-dimensional intersections is qualitatively the same as their behavior with (3 + l)-dimensional 

2 The authors of [IT] studied the D2-D8-D8 system at finite temperature but missed the existence of the long 
connected solution. 
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intersections. 





Figure 2: Behavior of £o//3 versus z* obtained numerically for generic values z* for various q's and intersections of 
interest, i.e. r = 1 and r = 3. The blue plots show the behavior for 3 + 1 intersections whereas the red plots show the 
behavior for 1 + 1 intersections. These plots show the number of connected Dp-branes (for all temperatures) placed 
in the background of black Dg-branes. For q = 5, the temperature T is fixed: T — (2tt.R5+i) _ . Except for q = 5,6, 
such solutions can potentially be realized as different phases of the holographic dual theories. 



4 Energy of the configurations 



As we saw in the previous section there are typically more than one solution for a given £q//3. 
In fact, just to recap, for £q//3 less than a critical value there are generically three branches of 
solutions for q < 4; one disjoint and two connected solutions. There are also an infinite number of 
solutions for q = 5 at T < as long as £q = 2irR^ + i/(r + 3). There are two branches 

of solutions for q = 6 for any £q/@; one disjoint and one connected solution. Since there are 
various configurations for a particular value oi£o/{3, one needs to compare their on-shell actions to 
determine which configuration is more energetically favorable. In this section we analyze the energy 
of these configurations by a combination of analytical and numerical techniques. The energy of 
these configurations by themselves is infinite. We regulate the energies of connected configurations 
by subtracting from them the energy of disjoint configurations 

e = umA^oo |y dzz*(i- -pfzr£) 2 - / dz zf / ' ( 26 ) 

where A is a cutoff and the difference in energy E is related to E as 

E = -T-C(q,r) u* / d r x E. (27) 
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D5 - D7 - D7 
at finite temperature 



D5 - D7 - D7 
at zero temperature 



1 <- z. 



z t — » oo 



<- </» 



-> oo 



D5 - D5 - D5 
at finite temperature 



D5 - D5 - D5 
at zero temperature 



<- y* 



Figure 3: The two plots on the left hand side show the number of connected solutions for the flavor D7 and 
D5-branes in the background of black D5-branes at T = (27r7?5+i) _1 . The plots on the right hand side, on the 
other hand, show the number of connected solutions for the same flavor branes in the background of thermal (or 
zero temperature) D5-branes. z*'s represent the turning points of the connected flavor branes in the black D5-brane 
geometry while y,'s are the turning points in the background of thermal D5-branes. 



4.1 Analytical analysis 

The integral in (|26p is complicated and as far as we know cannot be integrated analytically for 
generic values of z*. However, like the integral in (f2T|) . there are two regions of z* — > 1 and 3> 1 
for which we can integrate E(z*) analytically. In the z* — * 1 limit, it can be shown that E is 
positive for all q and r. Indeed, let's rewrite ([26]) as follows 




The last integral in ([28]) is positive for any z'* . If z* = 1 + e, we can choose z'^ = 1 + 3e and the 
difference in the bracket is estimated to be 

e(V6 -2) + e(log(\/2 + \/3) - l) +o(e 2 ), (29) 

which is also positive. 
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Figure 4: Behavior of E versus z* obtained numerically for generic values of z» for various q's. 



Approximating (|26p in the z* ^> 1 region yields 



E\ ~ J cte/ +7 (z 7 -z/) 2 - J dz 



jjT dz z 5+7 (V -zj) 2 -jf dzzf| + -^. (30) 



cj + 2 



Note that in the large z* limit, z* ~ zo where zq = uq/ut- As a result, from ([30]) we see that 
the energy of the connected configurations approaches their zero temperature value obtained in 
|12j . This is expected since in this regime the connected flavor branes are very far away from the 
horizon hence receive little effect from it. The disjoint configuration, on the other hand, always 
keeps in touch with the horizon and gets a finite temperature contribution —^2- The term in 
the curly bracket of (|30p is already computed in [12J in terms of Beta functions, giving an energy 
difference of 



E = limA^oo^A 

1 5 + 1 1 



±z} +S+1 B 

7 



2 7 ' 2 



+ ^T9- (31) 

a + 2 



Since z* S> 1, the last term is irrelevant in determining the sign of E. We see immediately that 
E > for q > 6 and E < for q < 4. For g = 5, the Beta function vanishes, and a more careful 
investigation must be made to determine the sign of E. It turns out that in going from (|26]) to 
(|30p . we have over-estimated the energy for the joined brane configurations: there is a correction 



with leading behavior ~ — z\ for large z*. Thus, for q = 5 we have £7 < for r > 1, while the 
analytic analysis is not reliable for r = 1. The analytic results for the two aforementioned limits of 
z* are summarized in Figure 4. 
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Figure 5: Behavior of E versus £o//3 obtained numerically for various q's and intersections of interest, i.e. r = 1 
and r = 3. The blue plots show the behavior for 3+1 intersections whereas the red plots show the behavior for 1 + 1 
intersections. The plots for the color D5-branes are for a fixed (inverse) temperature of /? = 2nR$+i. 

4.2 Numerical results and phase transitions 

For generic value of z* the integral in f)26|) can be numerically integrated and one can plot E versus 
z*. Instead we will numerically eliminate between ([21]) and ([26]) and plot E versus lo/(3. The 
reason for doing so is that one can easily observe the transitions, say from a connected configuration 
to a disjoint configuration, in terms of both t$ and (3. The numerical results are shown in Figure 
5. In the following we will analyze these graphs and determine the phases of the vacuum as one 
varies £o/(3- 

(3 + l)-dimensional intersections As mentioned in the previous section, for 3 + 1-dimensional 
intersections we have three configurations. Consider first the D4-D8-D8 system. For £o/(3 < 
(£o/P)t r ~ 0.154 the short solution is more energetically favorable compared to both the disjoint 
and the long solutions. This corresponds to the chiral symmetry being broken. (Here (^o//3)cr i s 
a critical value read off the energy diagrams in Figure 5. It is different from (£o/P)cr obtained 
from the plots of Figure 2.) Above this value (up to £q/(3 ~ 0.17) both short and long curved 
solutions have more energy compared to the disjoint solution, indicating that the disjoint solution 
is the vacuum, hence chiral symmetry is restored and the phase transition is of first order. Thus 
the chiral symmetry breaking-restoration phase transition occurs at (£q/P)^ t ~ 0.154. 

For the D5-D7-D7 system at T = (2it R^+i)' 1 the behavior is surprising and to some extent 
more involved. On the D5-D7-D7 plot in Figure 5 there are two special points with (^o)cr ~ 0.168 
and 0.17 (in units of 27ri? 5+1 ). For £ /(2irR 5+ i) < 0.168 the disjoint solution is more energetically 
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favorable, hence chiral symmetry in the dual theory is intact given that such a solution can represent 
a valid phase of the dual theory. For 0.168 < £q / '(2ir -R5+1) < 0.17 there are three kinds of solutions, 
disjoint, short and long. It turns out that the short solution has less energy than the other two 
indicating that chiral symmetry is broken. For £q/(2-kR5 + i) > 0.17 the disjoint solution becomes 
more energetically favorable, hence potentially chiral symmetry gets restored. For temperatures 
less than (27ri?5+i) _1 , the situation is the same as the zero temperature case: there are an infinite 
number of connected solutions when 3£o = 7ri?5 + i, each equally energetically favored, and each of 
them more favored over the disjoint solution. The existence of such solutions may be rooted in the 
fact that the low energy theory on the color D5-branes is a non-local field theory, a little string 
theory. Due to the fact that in the case of background D5-brane geometry the dual field theory 
degrees of freedom cannot be totally decoupled from non-field theoretic degrees of freedom, it is not 
clear to us whether such solutions can represent a chirally-broken phase of the dual gauge theory 
despite the fact that, geometrically, they are smoothly connected. 

For the D6-D6-D6 system, the disjoint solution is always favorable, although it is not clear 
whether one can give a holographic interpretation that the "dual" field theory is in a chirally- 
symmetric phase. This is because there exists no decoupling limit suitable for holography in the 
case of background D6-branes. 

(1 + l)-dimensional intersections For (l+l)-dimensional intersections, there are five allowed 
configurations. Among these the models based on background Dg-branes with q < 4 exhibit 
similar behaviors to their counterparts with (3+l)-dimensional intersection. That is to say for 
sufficiently low temperatures (compared to l/lo) chiral symmetry is broken while above a critical 
temperature it gets restored. The critical values at which this (first order) phase transition occurs 
are (£ /f3) c CT « 0.191, 0.196 and 0.206 for color D2, D3 and D4-branes, respectively. 

The model with color D5-branes again shows some surprises. Because of the mixing be- 
tween the field theoretic and non-field theoretic degrees of freedom in holography involving back- 
ground D5-branes, we have no evidence that different behaviors of the flavor branes represent, 
via holographic point of view, either chirally-symmetric or chirally-broken phases of the dual 
field theory. Nevertheless, one finds the following results. At T = (27ri?5 + i) _1 and for small 
enough £q (£0 < (£o)t T ~ 0.2498 x 27ri?5 + i) the disjoint solution is preferred. Increasing £q up to 
£0 1 '(2tt .R5+1) < 0.251 will result in a phase where the connected solution is favorable. This phase 
appears in our plot because we put a cutoff of A = 5 to regulate the energy integral. Increasing the 
cutoff will decrease the range of £q for which this phase exists. It is plausible that in the A — ► 00 
limit this phase disappears although our numerics does not allow us to check this explicitly. Sticking 
for now with the cutoff we chose, if one increases £q further, there will be another phase transition 
to a phase where the disjoint solution becomes favored. Due to space limitations, the resolution 
of the D5-D5-D5 plot in Figure 5 does not allow one to see all these phases. For temperatures 
less than (27nR5+i) -1 , like the zero temperature case, there are an infinite number of connected 
solutions for £q = 7ri?5 + i/2, each equally energetically favored, and none for other £q's. Each of 
these connected solutions is more favored over the disjoint solution. For the D6-D4-D4 system, 
there is no phase transition and it is always the disjoint solution which is energetically favorable. 
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5 Transverse intersections at finite temperature with compact x q 



An interesting property of the models we are studying here is that when x q is compact the scale 
of chiral symmetry breaking is generically different from the scale of confinement which results in 
additional phases. For example, for the Sakai-Sugimoto model at finite temperature, it was shown 
|10| that there exists an intermediate phase where the system is deconfined while chiral symmetry 
is broken. 



At finite temperature and x q direction being compact (with a radius of R c ), there are three 
geometries where the topology of the t — x q submanifold is S 1 x S 1 . One is a geometry which has 
the metric 



1=1 .._!=£., 



with 



9{u) = 1 -^my-\ (33) 

We will call this geometry the thermal geometry. Although the Euclidean time period (3 is arbitrary 
in this geometry, the x 9 -circle cannot have arbitrary periodicity. In order for this geometry to be 
smooth at u = ukk in the x q — u submanifold, one has to have 

7-q 

Az* = P C = ^- (^±L) 2 ukk, (34) 
7 - q V ukk ) 

where we have defined (3 C = 2nR c . Although we do not specify the (/-dependence of [5 C and R c , 
one should keep in mind that they depend on «kk and R q +i differently through (j34|) depending on 
what value for q is given. There is another geometry whose metric takes the form 

2 / „ „ „\ / II. \ 2 



^^{rTJ K + ^ + ^'j + ViT^J " { du2 + u2 ^s- q 2 ). (35) 

There is also the black brane geometry which is basically the same as ([2]) but with x q compact, 
and has the (Euclidean) metric 

* 2 = {it) V ( /(a)dt2 + * + + {it)^(m + A!V ' 2 ) ' (36) 

where 

/c>-i-(?r. or) 

In this geometry the a^-circle has arbitrary periodicity whereas (3 is fixed by 

fi=i£-(^j*«r. (38) 
7 — q V ut ' 

For all three geometries, the dilaton <f>, and the g-form RR-flux F q are given in ([6]); see the appendix 
for more details. Also, R q +i is given in (J4j). 
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The three geometries whose line elements are given in (|32p . (|35p and (|36p are saddle points of 
either type IIA or type IIB Euclidean path integral. For a given temperature, one needs to compare 
their (regularized) free energies to determine which solution dominates the path integral. We have 
calculated the free energies of these solutions in the appendix. For q < 4, both the thermal and the 
black brane geometries have less free energy compared to the geometry in (|35p . and this result is 
independent of temperature. So, to determine the lowest energy saddle point we need to compare 
the free energies of (j32|) and (|36p . The difference in their free energies is given by 

'-'thermal ~ ^black brane — ^ Vg \ f — q g+1 I { '> ■'• — "V ) • (•>>') 



s 



where Vg = /3/3 c Vol(S 8_9 )Vol(]R 9_1 ) where we set l s = 1. Thus the thermal geometry (I32p dominates 
when /3 > 27ri? c whereas for (3 < 2irR c it is the black brane geometry (|36p which gives the dominant 
contribution to the Euclidean path integral. This phase transition which happens at (3 = [3 C is the 
holographic dual of confinement-deconfinement phase transition in the corresponding dual gauge 
theories [6]. 

For q = 6, again, both the thermal and black brane geometries have less free energy compared 
to the geometry in (|35p . The difference in free energies for the geometries in (j32|) and (j36[) is given 
by 

^thermal D6 _ Sblack D6-brane = -. „ n ? (/3 2 — 0%), (40) 

where Vg = /3/3 c Vol(S 2 )Vol(M 5 ). So, ([30]) indicates a phase transition at (3 = 2ttR c . Unlike the 
q < 4 cases, for q = 6 the thermal geometry (|32|) dominates for (3 < 2ttR c whereas for (3 > 2irR c 
the black brane geometry (j36j) dominates. 



For q = 5, one needs to consider more possiblities. For /3 = 27ri?5 + i, there are two cases: either 
P c = 27ri?5 + i or f3 c / 27ri?5+i. For /5 C = 27ri?5 + i, both ([32]) and (|36|) are more favored over the 
geometry whose metric is given in (|35p . The difference in free energies of (|32|) and (|36|) is 

^thermal D5 _ Sblack D5-brane = ^^(it 2 -. — n|; K ), at f3 c = 2irR 5+ i. (41) 

9 s 

where Vg = /3/3 c Vol(S 3 )Vol(E 4 ). For (3 C ^ 27ri? 5+1 , on the other hand, the two saddle points with 
the same asymptotics are the black brane geometry and the geometry in (|35p . In this case, the 
black brane geometry is the dominant one 

AS = ^V 9 ul > 0. (42) 

There are also two possibilities when (3 / 2t^R^ + \. If (3 C ^ 27ri?5 + i, the only saddle point consistent 
with the asymptotics is (|35|) which determines the vacuum. If, on the other hand, (3 C = ZnRs+i, 
there are two geometries with the same asymptotics, thermal and the one given in (|35p . The 
thermal geometry is dominant because 

AS = -^V 9 u 2 KK < 0. (43) 
9% 



With x q being compact, the flavor branes are now sitting at two points separated by a distance 
£q on the x^-circle. Note that there is no reason for the flavors branes to be located at the antipodal 
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points on the circle. As before, we consider the transverse intersections of the flavor and the color 
branes, and choose the same embeddings and boundary conditions for the flavor branes as we did 
in (fl"0|) and (fTTj) except that now £q < ttR c . In what follows, we will focus on q < 4 cases and 
consider both low and high temperature phases of the background. In particular, we would like 
to know whether there always exists a range of temperature above the deconfinement temperature 
where chiral symmetry is broken. 

5.1 Behavior at high and low temperatures 

For temperatures above the deconfinement temperature (3 C the profile of the flavor branes takes 
essentially the same form as it did when x q was non-compact, hence, indicating the existence of 
short and long (smoothly) connected solutions. Note that when x q is compact, for a fixed 1$, there 
is now a lower bound on £o/P set by the deconfinement temperature. For completeness, we have 
plotted £o/P versus (z* being the radial position at which the brane and anti-brane smoothly 
join) in Figure 6. The lower dotted line in each plot represents the deconfinement temperature. 
As an example, we chose it to be at (3 C = lO^o- The upper dotted line shows chiral symmetry 
breaking-restoration phase transition which comes from comparing the energies of connected and 
disjoint solutions. One can show, using energy considerations, that below the upper dotted line 
chiral symmetry is broken in a deconfined phase via short connected solution while it is restored 
above the line (where we have deconfinement with chiral symmetry restoration). 

D4-D8-D8 D2-D8-D8 




Figure 6: Behavior of lo/P versus z, above the deconfinement temperature for q = 2,3,4. The lower dotted line 
represents a transition to the confined phase whereas the upper dotted line represents chiral symmetry breaking- 
restoration phase transition. 
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In the low temperature regime where the system is in a confined phase, the DBI action for the 
flavor branes in the thermal background (|32p now reads (with gauge fields set equal to zero) 



'DBI 



(3C(q,r) / d r x dx q U2 g(u) + 



1 



U \ 2<5 /2 
U 



(44) 



g(u) \Rg + i- 

where C(q, r), 7 and 5 have all been defined in (|14p . The equation of motion for the profile is now 

1 



U-2g(u) 



g{u) + 



g(u) \R, 



U \ 2S i2 

) U 

•9+1 ' 



-I ^ 

Z 9 

= w. 



- 



(45) 



with wo parameterizing the solutions. There exists a solution with wo = representing disjoint Dp 
and Dp-branes descending down to u = ukk- For wq ^ 0, solving (l4"5j) for it yields 



1 



-25 / \ 



(46) 



Denoting the possible turning point (s) by w*, analysis of (f4~6"j) shows that u > w* > ukk, with w* 
satisfying 



where a has been defined as before. Integrating (|46p gives 



-25 



1 



(47) 



(48) 



where we have defined y = (u/ukk) G (1>°°)> an d 2/* = ^*/^kk- Using (|48l) we can relate y* to £q 

-1 



^ = — (yl-y*)^ 

Pc K 



y- 25 - 1 



y 1 -y a -{y2-yl)) dy. 



(49) 



The analysis of £0/ ' f3 c as a function of determines the number of solutions. £0/ (3 C versus has 
been numerically plotted in Figure 7 for intersections of interest and for q = 2, 3, 4. As it is seen 
from Figure 7, there is always one smoothly connected solution (as well as a disjoint solution). 
The red and blue plots represent smoothly connected solutions for (l+l)-dimensional and (3+1)- 
dimensional intersections, respectively. 



For the energy of the connected solutions, one obtains 

E = -^-PMq, r) 4k / dx r (50) 

where 

E = hm A _ { [ A dy yi (l - y^) ~ k (l - ^4) " t dy (l - y") " V } , (51) 

and A is a cutoff. Like the previous sections, one can numerically eliminate y* between ()49|) and (|5ip 
and plot E versus £o/Pc- Although we have not shown the plots here, one can check (numerically) 
that the smoothly connected solutions are always more energetically favorable compared to the 
disjoint solutions. 
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D4-D8-D8 D2-D8-D8 




Figure 7: Behavior of £ //3 c versus w, for q = 2, 3,4. 



The plots in Figure 7 not only show the existence of a unique smoothly coonected solution 
(for the confined phase) but also indicate a big difference for the behavior of the flavor branes 
below and above the confinement-deconfinement phase transition. For example consider the plot 
for D3- D7-D7 system. Each point on the plot represents a unique curved solution for which 
increasing the temperature (up to the deconfinement temperature /3 C ) will have no effect on the 
shape of the U-shaped flavor branes. That is to say that changing the temperature will not cause 
the flavor D7-D7-branes to join either closer to ukk or farther away from it. Once £o and (3 C are 
specified, the shape of the brane stays the same independent of the variation of the temperature up 
to the deconfinement temperature. The behavior just mentioned is significantly different from the 
behavior of the flavor branes at high temperatures (where they are in black brane backgrounds). 
For a fixed £q, each point on the plot of D3-D7-D7 system in Figure 7 represents one (or two) 
smoothly connected solution(s) for only a specific temperature. Varying the temperature will now 
change the shape of the connected flavor branes and force them to go either closer to the horizon 
or stay farther away from it. 

There is another difference which is worthy of mentioning here. At the temperature for which 
the confinement-deconfinement phase transition occurs, namely at (5 = /3 C , ukk = ut- For a fixed 
lo, flavor branes at temperatures above /3 C join at a radial point closer to ukk = ut than the same 
flavor branes placed at temperatures below f3 c . 
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6 Discussion 



In this paper, we analyzed some aspects of transversely-intersecting Dg-Dp-Dp-branes at finite 
temperature. In particular, we mapped out different vacuum configurations which can holographi- 
cally be identified with chiral symmetry breaking (or restoration) phase of their holographic dual 
theories. Although we showed that generically the long connected solutions are less energetically 
favorable compared to the short connected solutions we did not discuss their stability against small 
perturbations. Presumably a stability analysis along the lines of [19j can be done to show that the 
long connected solution is unstable against small perturbations. The analysis presented here can be 
generalized in various directions. For example, one can add a chemical potential to the setup and 
look for new phases as was done in some specific models in [20l EU [22J, [23] , or consider the sysytem 
at background electric and magnetic fields |24[ |2"5[ [26] and study the conductivity of the system or 
the effect of the magnetic field on the chiral symmetry-restoration temperature. We hope to come 
back to these interesting issues in future. 

Also, our analysis was entirely based on the DBI action for the flavor Dp-Dp-branes. In 
transversely-intersecting D-branes, working with just the DBI (plus the Chern-Simons part of 
the) action misses, from holographic perspectives, an important part of the physics, namely the 
vev of the fermion bilinear as an order parameter for chiral symmetry breaking. To determine the 
fermion bilinear from the holographic point of view there should be a mode propagating in the 
bulk geometry such that asymptotically its normalizable mode can be identified with the vev of the 
fermion bilinear. The DBI plus the Chern-Simons action cannot give rise to the mass of the localized 
chiral fermions either. Note that for transverse intersections one cannot write an explicit mass term 
for the fermions of the intersections because there is no transverse space common to both the color 
and the flavor branes. So it is not possible to stretch an open string between the color and flavor 
branes in the transverse directions. Therefore, the fermion mass must be generated dynamically. 
It has been argued in [8] that including the dynamics of an open string stretched between the 
flavor branes into the analysis will address the question of how one can compute fermion mass and 
bilinear vev in these holographic models. More concretely, the scalar mode of this open string which 
transforms as bifundamental of U(N{) x U(N{) has the right quantum numbers to be potentially 
holographically dual to the fermion mass and condensation. 

Recently, the authors of [28j [29j [30] shed light on this issue by starting with the so-called 
tachyon-DBI action [27] claimed to correctly incorporate the role of the open string scalar mode, 
the tachyon, in a system of separated Dp-Dp-branes. In fact, it was shown in [29, 30J that for 
the Sakai-Sugimoto model at zero temperature, the open string tachyon will asymptotically have a 
normalizable as well as a non-normalizable mode. They identified the normalizable mode with the 
vev of the fermion bilinear (order parameter for chiral symmetry breaking) and the non-normalizable 
mode with the fermion mass. It is not hard to generalize the calculations of [291 [30] to include all 
transversely-intersecting Dg-Dp-Dp systems at zero temperature where one finds that there always 
exist both normalizable and non-normalizable modes for the asymptotic behavior of the tachyon, 
and in the bulk of the geometry the tachyon condenses roughly at the same radial point where 
the flavor branes smoothly join |31j. The calculation of the fermion mass and condensate for 
transversely-intersecting Dq-Dp-Dp systems at finite temperature requires not only considering the 
tachyon-DBI action of [27] in the black brane background of ([2]) but also calculating the tachyon 
potential as a function of the temperature. In the case of a coincident brane and anti-brane in flat 
background, a partial result for such a calculation was given in [32] (see also [33]). Some attempts in 
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generalizing the results of \32\ [33] for separated brane-anti-branes (at least in the case of separated 
D8-D8 in flat space) has recently started in [33]. For our purpose of extracting information about 
fermion mass and condensate, knowing the dependence of the tachyon potential for separated flavor 
Dp-Dp-branes seems crucial. Ignoring the temperature dependence on tachyon potential at zeroth 
order yields unsatisfactory results: It gives rise to the same results as one would have obtained 
in the zero temperature case |31j . We know that this is not the right behavior because at zero 
temperature there is no chiral phase transition. 



Note added in version two There are now two more interesting (and closely related) proposals 
on how to compute the fermion mass and condensates in transversely-intersecting D-branes. One 
proposal is based on an open wilson line operator whose vev gives the condensate [35] (see also 
[361 [37] for generalizations) and the other [38] is based on a D6-brane, in the case of the Sakai- 
Sugimoto model, ending on the flavor D8-D8-branes It would be interesting to understand the 
relations between these two proposals and the one based on the open string tachyon. 
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A Free energies of possible bulk backgrounds 

Following [10] we use the notation adopted in [39] to compute the free energy of possible finite 
temperature Dq-brane backgrounds for both compact and non-compact x 9 -direction. We first 
consider the case of compact x q and as a warm-up first do the computation for D6-branes, then 
generalize it to Dq-branes. 



A.l Compact x q 

D6-branes Consider the bosonic part of type IIA supergravity action (in Euclidean signature 
and) in string frame 

S = — (Seh + Stj, + Srr + 5ns), 

= -J e- 2 ^(TZ + 4W) + ~ Y, j F (i+2) A * F (i+2) + \J e ~ 2< ^(3)- A *H {3) . (52) 
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Having the near horizon geometry of D6-branes in mind, we choose the following (radial) ansatz 
for the metric and RR 7-form C(j\ 

l~ 2 ds 2 =dr 2 + e 2X ^dx\ + e 2X ^dx 2 c + e 2 ^dn 2 , (53) 
C (7 ) =A dx 1 A ... A dx 5 A dx 6 ~ c A dx c , (54) 

where dVL 2 is the line element of the unit two-sphere. Certainly the backgrounds we considered in 
this paper can all be put into the form of our ansatz fj53[) . Here r is the radial coordinate, xn are 6 
(out of 7) worldvolume directions, x c plays the role of time in black D6-brane geometry and x 6 in 
thermal geometry. From now on we set l s = 1. It could be restored into our results by dimensional 
analysis. 



The equation of motion of C(y\ has the solution 



where the dot represents the derivative with respect to r. Substituting (|55j) into (152]) . the RR part 
of the action reads 



-S RR = J Q 2 e M+ ^dr, (56) 

where Q 2 is an integration constant. It is more convenient to define a new radial coordinate 
dp = —e^dr, where 

ip = 2<p - 6A - A - 2u, (57) 

in terms of which (|56j) reads 

-S RR = -J Q 2 e 6X+ ' x - 2v ^dp. (58) 
Note that ■ s /ge~ 2l ^dT = —e~ 2ip dp by which the dilaton action is easily computed to give 

- = f A<p' 2 dp, (59) 



where prime denotes the derivative with respect to p . With our metric ansatz (J53J), the Ricci scalar 
is 

n = 2e~ 2u - 42A 2 - 2A 2 - 6z> 2 - 12AA - 24Az> - 4Az> - 12A - 2A - 4i>. (60) 
Thus, expressed in terms of the p coordinate, the Einstein-Hilbert action becomes 

- S'eh = J Ke~ 2ip dp, 

= j \2e~ 2u ~ 2v - 42A' 2 - 2A' 2 - 6i/ 2 - 12A'A' - 24AV - 4AV - 12(A" + \'<p') 

-2(A" + XV) - A{u" + v'ip')\ dp. (61) 

Note that the terms with double primes are total derivatives and are cancelled by adding the 
Gibbons-Hawking term to (f6T|) , Using ([57]) . S'eh and S^ add up to give the following simple 
expression 

- S EH - S^ = J (2e~ 2v ~ 2 * - 6A' 2 - A' 2 - 2^' 2 + c^' 2 ) dp. (62) 
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Putting everything together, the total action reads 

S = f ( 2e~ 2 ^ - 6A /2 - A' 2 - 2z/ 2 + p' 2 - Q 2 e 6X+ ~ x - 2 ») dp. (63) 



The equations of motion are 



A" = ^! e -2*+12A+2A ) (64) 



A" = ^! e -2*+12A+2A ) (g5) 



J, 3Q 2 _ 



20+12A+2A 



(67) 



2 

where ip has been replaced by 4> using ([57]). Defining = — 2</> + 12A + 2A, ([64]) . (|6"5j) and (|6T|) give 

$" = 4Q 2 e*, (68) 

with the following solution 



$ = -2 In (^Q- sinh C lP ) - 2 In C. 



(69) 



The last term in (|69p is there just for convenience. Going back to (|64p. (|65p and (|67p and solving 
for A, A and 4>, one obtains 

A = - \ In smhC7 1/0 ) + C 2 A p, (70) 



4 V 5s Ci 

A = - ~ In sinhC 1/0 ) + C x p, (71) 

^ = -7ln(^sinhCip) +C|p + ln 5s , (72) 

with 6C2 + C*2 — = 0. The constants C±, C^, and are to be determined. Defining 
$0 = 3> + 2 In the equation of motion for ^ becomes 

y" = e ^-2C$p-41n gs+ 2v _ l g 2 e * (?3) 

If we take 2v = + 2C2 p + C^, the above equation is satisfied if 

C u = \n2Q 2 g 2 s . 

Thus, we obtain 

v = A\ n sinhCip) + C+p + lnV2Qg s . (74) 

4 Vg s Ci / 
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In what follows we show that the near horizon geometry of non-extremal D6-branes satisfies 
(HH, CD}, ((72]). This geometry takes the form 

ds 2 =f~ 1/2 [dx\ + hdxfj + f 1/2 {^du 2 + u 2 dn 2 2 , ) (75) 



-2c£ __-2 

with 



e ^ = 5s V 3/2 , (76) 



/ = |^, and h=l-^, (77) 

where u c is either or -ukk (or zero for the geometry given in (|35p ). The u coordinate is not the 
same as p. Their relation can be inferred from (|70p and (jTlf) . and is given by 



yielding 



2u 2 AC 2 h 

where AC2 = C 2 — C 2 . Computing e 2A from (|70p . we have 



= 2(A-A) = e 2p(C}-C-) = K (7g) 



d P = o-p"a^ T dn > ( 79 ) 



e 2A 



@)sinhC lP 



c. 



/j2AC 2 _ /j2AC 2 



C 2 __£l 



<3 / E^T 
1 - h A °2 



(80) 



To match e 2A in ([80]) with Z 1 / 2 , one has to impose —C\ = AC2 and Ci = 4C*2 , or 

Ci = 4C 2 A ; C\ = -3C 2 A. (81) 

Thus, we have 

„2A _ / V2C\g s ~u — ( 2l 3 u \ 2 . _ Qu c l s 

e - V O = UjV =^ °2 - n,/5„2W 



v 9^y 2 2v / 2 ff fjV 

Taking = 36*2 , we can match the dilaton with the dilaton of the D6-brane solution (|76|) 

e 2, = e 6A + 21n 9s=9 2/^ 3 / 2 _ (g3) 
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Furthermore, from (|74p . we have 

° 2 "=^=^4M <2 -&" 2 ^ N = 2 ^ Ql ° (84) 

The final check is the g uu component of the metric. To do that, we first need to compute e~ 2ip 
which turns out to be 

e -2tp _ g -40+12A+2A+6i^ 



e 



14A-8C*^p+4 In ( V2Qg a ) -4 lng s 



2lsU ) V \(V2Q)\ (85) 



Hence, 



dr 2 = e-^dr 2 



P 



7 4u 4 (l6(C 2 A ) 2 J 71 



■v / - ( 86 ) 
2i s u / h 

where we have used (I82p and (1841) to simplify the expression. Thus, the near horizon geometry of 
N c D6-branes is indeed a solution to type IIA supergravity equations of motion. 

It is now straightforward to plug ([70]) . ([71]) . (172]) and ([73]) into ([63]) to obtain the on-shell action 
(free energy). Changing the coordinate from p coordinate to u and being careful about the fact 
that dp/du is negative, we obtain 

S = _y 9 jT {(- 6A' 2 - A 72 -2i/ 2 + V ' 2 )^+^^ (87) 

= \v 9 / A du = 4 y 9(A - «c). (88) 
9s J u c 9s 

where in the first line prime denotes the u-derivative, and Vg is the volume (in units of string length 
which we set equal to one) of space transverse to u: Vg = /JVol(Si q )Vol(S 2 )Vol(K 5 ). The free energy 
is infinite. That is why we put a cutoff u = A in the second line to regulate the free energy. The 
cutoff drops out when it comes to comparing the free energies of different solutions. 



It is clear from (|87p that a D6-brane solution of the type (|35h has always more free energy 
than a black D6-brane (|36]) or a thermal D6-brane (|32|) . Comparing the free energies of the black 
and thermal D6-brane geometries, we obtain 

3 

^thermal D6 — Sblack D6-brane = —^Vg(uT ~ Ukk) 

9% 

3 ^9 I n2 ,1 



l6TT 2 g 2 RQ +1 

where in the second line we have used ([33]) and ([38]) . and 

47T 



(P-K), (90) 



C = 2ttR c = — v/i? 6+1 UKK- (91) 

Thus, there is a phase transition at (3 = 2ttR c , where the thermal D6-brane geometry dominates 
for (3 < 2-7rR c whereas for f3 > 2nR c the black brane geometry dominates. 
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Dg-branes: general analysis Having gone through the analysis for the D6-brane, calculating 
the free energy of non-extremal Dq-branes is straightforward. The relevant parts of the supergravity 
action (in string frame) are 

S = — (5eh + S<t> + Srr + <Sns), 

= - j e- 2 ^(K + AdcpdcP) + ^ E / F (i+V A *F (q+2) + i J e~ 2 *H {3) A *H {3) . (92) 

We are interested in backgrounds with Hrg\ not turned on, and choose the following ansatz for the 
metric and the RR g-form C g +i 

ds 2 = dr 2 + e 2X ^dx\ + e 2 ' x ^dxl + e 2u ^dnl (93) 
C {q+1) = A dx 1 A ... A dx q ~ x A dx q ~ c A dx c , (94) 

where k = 8 — q, and dQ 2 . is the line element of unit fc-sphere. Again, we have set l s = 1. The 
equation of motion for C( g +i) has the solution 

Ae -gX-X+ku = congt (95) 

The RR action then becomes 

- S RR = [ Q 2 e" x+ ~ x -^dT, (96) 



where, again, we denoted the integration constant by Q 2 . For convenience, we change the variable 
from t to p defined by dp = —e^dr, where 

tp = 2<t> - qX - A - kv. (97) 

Then 

- S RR = - [ Q 2 e qX+x - k ^dp. (98) 



Following the same argument as we did for q = 6, it is not hard to see that the action (|92p 
takes the form 

S = [ (k(k- l)e~ 2v - 2 ^ - q\' 2 - X' 2 - kv 12 + ip' 2 - Q 2 e qX+x ~ ku ) dp. (99) 



Having put the action in the above form, the equations of motion are now in order 

A" = ^! e -20+2 g A+2A (1Q0) 

A" = ^! e -20+2 9 A+2A (1Q1) 

v » = ( k - l)e- 2 "- 2 * - ^! e -2^+2 g A+2A (1Q2) 

^ = (5 - k)Q 2 e _ 2 ^ +2gA+2 A ; (103) 
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Define $ = -20 + 2qX + 2A, (fTOOl) . then (fTOTT) and (fT03l) result in 

= 4Q 2 e*, (104) 

which yields the following solution 

$ = _21n(^^sinhCip) -21nC. (105) 

The last term is there just for convenience. Solving for A, A and (p in (jlOOp . (jlOip and (|103p . we 
obtain 

A = Ai Q (^Q- sinhCl p) +C£p, (106) 

A = _Ii n (^ sinh C 1/0 ) +C 2 V (107) 

^ = -^ln(^sinhC 1/3 )+^p + ln 5s , (108) 

with qC$ + C| - C| = 0. The constants Ci, C£, c| and Cf are to be determined. Define 
$o = $ + 21n<7 s , then the equation of motion for i/ becomes 

^ = (fc _ l)e ^*o-2C^-41n ffs+ 2(fe-l), _ ^_ Q 2 e * (1Q9) 

If we take 2(k — l)i> = ^^^o + 2C^p + C„, then the above equation is satisfied if 



Thus, we obtain 



5-fc /V2Q . n Cfp 1 ^Qff, 
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We now check whether the near horizon geometry of non-extremal Dq-branes satisfies (|106p . 
(fT07l) . (fT08l) and (fTTTI) . This geometry takes the form 

ds 2 = f~ 1/2 (dx\ + /idx 2 ) + fl 12 {^du 2 + u 2 dn 2 k ^j , (112) 

e- 2 * = <77 2 4 9 " 3)/2 , (113) 

where f g = ^ q9 ^- q , d' p = (^V^) 5 ~ 9 ^('^-) an d fa = 1 — {u c /u) 7 ~ q .The relation between u and pfao 
can is inferred from (|106p and ()107p 

e 2(A-A) = e 2p(c|-^) = ^ (n4) 

which gives 
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where AC2 = C 2 — C 2 . Calculating e 2A from (|106p . we have 



e 2A 



e 2C 2 V 



(V2Q\ 

\asC\ 



s'mhCip 



c. 



! V2C l9s 



Q 



h2AC 2 _ / l 2AC 2 

C 2 _ Oy 



" V v / ^ ( } 

The only way to make e 2A in (|116p to be equal to is by imposing — C\ = AC2 and Ci = 4C 2 , 
or 

Ci = 4C 2 A ; C 2 A = -3C 2 A. (117) 

Thus, we have 



e 



2X= V2c l9s ,u^^ 1/2 A = qy s ) 7 -« (118) 

\ Q \u c J Jq 2 4V2glN c d' g K ' 



J2<t> _ p 2A(5-fc)+21n 9s = n 2 f- 3 ^ 



Also, the dilaton can be matched to (jll3[) as long as we take C 2 = (5 — fc)C 2 resulting in 

9lU^. (119) 
Furthermore, from (lllip . we have 

e^ = (^)^ T e 2A ^^/ g 1 /V => N c = (120) 
The final check is the g uu component. To do that, we first need to compute e~ 2ip which gives 

-2<p e -i<t>+2q\+2\+2kv 

= e ^-^P+^ln(^)~4ln 9s _ (121) 



Hence 



dr 2 = e^dr 2 



P 

{7-q) 2 u 2 c (7 - q) u 2 ^) 1 6M^ A _ 8C A p+#rln( y|Q, 



64(C A ) 2 h 2 



= / g 1/2 ^^ 2 , (122) 

where we have used (|118p and (|120p to simplify the expression. Therefore, the exact solution we 
just found indeed corresponds to the near horizon geometry of non-extremal N c Dg-branes. 
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With ([T06]) . (fTUT]) . (fT08|) and (fTTT|) at hand, it is now straightforward to plug them into ([99]) to 
get the on-shell action. Changing back to the u coordinate and being careful about the fact that 
dp/du is negative, we obtain (here prime denotes u derivative.) 



S = -y 9 |J { ( - 9 A' 2 - A' 2 - b' 2 + V ' 2 ) | + - l)e" 2v -^ - Q 2 ^"" j ^' 

'-Vg I (7-q) U 6 ~ q du 



■du 



9-g. 



9 2 s 

= ^S-V 9 (A 7 -i-u 7 -i), (123) 

where A is the cutoff, and Vg = (3R c Vol(S 8 ~ q )Vol(M. q ~ 1 ) (in units where string length l s = 1). For 
q < 4, it is not hard to see from (|123p that the geometry given in (|35p for which u c = has more 
free energy than the geometries given in (|32p and (|36p . The difference in free energies of the thermal 
and black Dg-brane geometries is given by 

^thermal ~ Sblack brane = ^ 2 ~ ^9^T ? ~~ %k)' (124) 

which using (f34j) and (f38j) can be equivalently expresses as 

Q Q 9-g / 4^ I( 7 _ 9 )\ 2 H / 2 tl 2§^\ 

'-'thermal ~ ^black brane — ^ 2 ^9 I ^ _ ^ g+1 J I " _ " c J ' l i2t) J 



Thus, for q < 4 there is a phase transition at j3 = 2ttR c . For /5 > 2irR c the thermal geometry (|32p 
dominates whereas for /3 < 27ri? c it is the black brane geometry (j36[) which dominates the Euclidean 
path integral. The case of D5-branes and any associated phase transition has been discussed in 
section 5. 



A. 2 Non-compact x q 



When x q is not compact, and g / 5, there are only two geometries with the same asymptotic 
boundary condition: the thermal geometry ([7]) and the black brane geometry ([2]) . The difference 
in free energies of the two geometries is obtained by first setting u c = in (|123l) for the thermal 
and u c = ut for the black brane geometry, then subtracting the results 

^thermal — S'black brane = J^VgvL 9 , (126) 

9% 

where Vg = Vol(S 8_9 )Vol(IR' ? )Vol(S^) in units where l s = 1. The difference in free energies ([8]) is 
positive indicating that independent of temperature the black brane geometry is always dominant 
and determines the vacuum. The case of D5-branes has been discussed in section 2. 
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